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Abstract 

The Euler-Maxwell system describes the evolution of a plasma when the collisions are im- 
^^ ' portant enough that each species is in a hydrodynamic equilibrium. In this paper we prove 

Cn I global existence of small solutions to this system set in the whole three-dimensional space, by 

combining the space-time resonance method, dispersive estimates, localization estimates and 
energy estimates. An important novelty is that we can prove a very slow growth (t'^'^) of high 
derivatives even with only -K= decay of the L°° norm by reiterating the energy estimate. 
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1 Introduction 

1.1 Plasma physics and Euler-Maxwell 

There are different models to describe ttie state of a plasma depending on several parameters such 
as the Debey length, the plasma frequency, the collision frequencies between the different species... 
Formal derivation of these models can be found in Plasma Physics textbooks (see for instance 
Bellan [1], Boyd and Sanderson [3], Dendy [8] and the paper [2] ...) 

Since the plasma consists of a very large number of interacting particles, it is appropriate to 
adopt a statistical approach to describe it. In the kinetic description, it is only necessary to evolve 
the distribution function fa{t,x,v) for each species in the system. Vlasov equation is used in this 
case with the Lorentz force term and a collision term. It is coupled with the Maxwell equations for 
the electromagnetic fields. 



If collisions are important, then each species is in a local equilibrium and the plasma is treated 
as a fluid. More precisely it is treated as a mixture of two or more interacting fluids. This is the 
two-fluid model or the so-called Euler-Maxwell system. We refer to [23t \T6[ [2T] for more about 
hydrodynamic limits. Another level of approximation consists in treating the plasma as a single 
fluid by using the fact that the mass of the electrons is much smaller than the mass of the ions. 
This is the model which we are going to consider in this paper. 

1.2 The Euler-Maxwell equation 

The Cauchy problem for the one fluid version of the Euler-Maxwell system reads 
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The unknown functions are: /), the density of electrons; u, the average velocity of the electrons; 
E, the electric field; B the magnetic field. The physical constants are: c, the speed of light; e, the 
charge of the electron; tti, the mass of the electron. Finally, p is the uniform density of ions, and 
the electron gas is supposed to be barotropic, the pressure being given by p{p). 

Let us first recall a few results related to (jl.ip . Global existence of weak solutions was obtained 
for a related Id model in [5] using compensated compactness. Also, several asymptotic problems 
(WKB asymptotics, incompressible limit, non-relativistic limit, quasi-neutral limit...) were studied 
to derive simplified models starting from the Euler-Maxwell system |3H [33l W2\ [28] . We also refer 
to |25] where the incompressible Navier-Stokes system is studied. 

Going back to our system (jl.ip . we notice that the two last equations above can be removed, 
as soon as they are satisfied at the initial time, which we assume from now on: they are then 
conserved by the flow given by the first four. 



1.3 Vicinity of the trivial equilibrium state 

An obvious equilibrium state of the above system is {p,u,E,B) = (/>, 0,0,0). In order to study 
its stability, it is instructive to linearize the above system, and compute evolution equations for its 
unknowns. It is convenient then to split u into its divergence- free part Pu, and its curl- free part 
Qu: u = Pu + Qu. Similarly, E = PE + QE. The obtained system reads 



{d! - clA 



, dt{B- ^V X u 
where the speed of sound Cs and the plasma frequency ujp are given by 
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Thus around the equihbriuni, and at a hnear level, some unknowns are governed by Klein Gordon 
equation (with different speeds), whereas the quantity B— ^^V xuis conserved. The Klein Gordon 
equations entail decay, which is one of the keys of the global stability result which we will prove; 
as for the quantity B — ^^V x u, no decay is to be expected a priori. We will therefore set it to 
zero, which, as it turns out, is conserved by the nonlinear flow. 

1.4 Adimensionalization and reductions 

In the following, we set for simplicity the physical constants m, e, c, as well as p to 1. We also drop 
the 47r factors, since they are irrelevant. However c^ = p'{p) = ^'(1) remains a number less than 1. 
In order to simplify a little bit the estimates, we assume 



Finally, set 

The Cauchy problem becomes 



pip) =' |p^. 



def ^ 
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dtu + u ■ \7u = —cIpVp — E — u X B 

dtp + V ■ {pu) = 

dtB + V X E = {) 

dtE-V xB = pu 

\/ . E = -n 

V ■B = 

{u,n,E,B){t = 0) = {uo,no,Eo,Bo). 



We shall furthermore assume that, initially. 



B = V xu. 
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This condition is conserved by the flow of the above system: in order to see this, use the identity 
u ■ Vu = —u X (V X n) + V^ to compute 



dt{B - V X u) = V X {u ■ Vu + u X B) 

= V X (-ux{Vxu)+V- 



V X {ux B) 



The linearized system reads now 
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1.5 Obtained results 

Prior to stating our theorem, we need to define the operator A = |^ (see Section [2] for the precise 
definition of this operator). 

Theorem 1.1. Assume that the resonance separation condition \4-l\ holds; it is the case generically 
in Cs- Fix qq > 0. Then there exists Co,eo,NQ > such that: if e < eo, N > Nq and 

then there exists a unique global solution of (EM) such that 



sup 

t 



{t)-'''^'\\{u,An,E,AB){t)\\HN + y^\\iu,An,E,AB)m3 



< 



Furthermore, it scatters as t goes to infinity in that there exists a solution (ui,ni,Ei,Bi) of the 
linear system lll.&\) corresponding to intial data in H^^^ such that 

\\{u,n,E,B){t) - {ui,ni,Ei,Bc){t)\\jjN-2 -^ as t^ oo. 

Remark 1.2. A few observations on the hypotheses on the initial data: 

• The requirements on Auq and ABq imply necessarily that f n^ = and f Bq = 0. In particular 
this is consistent with the electric neutrality. 

• We did not try to optimize the number of derivatives in L^ required (N), but rather aimed at 
a proof as simple as possible. On the other hand, the weight appearing above ({x)^~^°) seems 
nearly optimal; a more precise analysis would maybe allow (x) instead o/(x)^^". 

The proof will be essentially split into two parts: controling the H^ norm of {u,n,E,B); and 
proving the decay in various norms. The former is achieved by an energy estimate; and the latter 
by the method of space-time resonances, which was introduced in |12j . It was also used to prove 
global existence of small data solutions for water waves [13t [H] . 

1.6 Stability of compressible Euler and related models in dimension 3 

It is instructive to compare the above results to earlier works on compressible Euler in dimension 
3, or couplings of compressible Euler with various fields (electrostatic, electromagnetic, gravita- 
tional...). For all these models, a fundamental question is whether given data lead to blow up or a 
global solution. 

A first class of results gives blow up for various types of data. The fundamental work is due to 
Sideris [30j, who proved finite time blow up for compressible Euler; he was able to obtain this result 
for data arbitrarily close to the equilibrium state given by a zero velocity, and a constant density. 
Many results followed: finite time blow up was showed for the compressible Euler equation with 
compactly supported data by Makino, Ukai, and Kawashima [22]; for the attractive Euler-Poisson 
equation with compactly supported data by Perthame [27J; for the repulsive Euler-Poisson equation 
with compactly supported data by Makino and Perthame [23] ; and for the relativistic compressible 
Euler equation by Guo and Tahvildar-Zadeh [19] and Pan and Smoller [26]. 

All of the aforementioned results rely on a non-constructive proof, and do not say much about 
the nature of the singularity. Recently, Christodoulou [6j was able to describe in a very precise 
manner the blow up process for the relativistic compressible Euler equation. 



Another line of research gives global existence (and scattering) for data close to the equilibrium 
state given by constant density, and all the fields (including the velocity) equal to zero. Such of 
result was first obtained by Guo |17j for repulsive Euler-Poisson; and by Guo and Pausader |18J for 
the ion dynamics in Euler-Maxwell. In both cases, the curl of the data is assumed to be zero, and 
this condition is conserved by the flow of the equation. Finally, global existence for Euler-Maxwell 
with relaxation was obtained by Duan [9]. 

Focusing on the case of small data (i.e. close to an equilibrium), some common features emerge 
from the results which have been mentioned. Global existence is only known under the assumption 
that the flow is irrotational: this eliminates a mode which is linearly non-decaying. Under this 
assumption, a crucial point is then the nature of the linearized equation: roughly speaking, blow 
up may occur if it is a wave equation, whereas global existence is expected if it is a Klein-Gordon 
equation. The relevant difference between these two situations is that the latter gives a decay 
~ -i, whereas the former only decays ~ j. 

In the case of Euler-Maxwell, which is treated in this paper, the condition i? = V x n is also 
meant to restrict the solution to the subspace along which the linearized problem is governed by 
Klein-Gordon equations. The novelty is that these Klein-Gordon equations have different speeds, 
making the nonlinear interaction more intricate. 

2 Notations 

We shall use the following standard notations: 

• A < B ii A < CB for some implicit constant C. The value of C may change from line to line. 

• A^ B means that both A<B and B < A. 

• For any real number a, the "Japanese brackets" (•)q, stand for {x)a = \/l + a^x"^. We also 
denote (x) = (x)i. 

• If / is a function over M^ then its Fourier transform, denoted by /, or Ff, is given by 

m = ^m = (^ / -~'''f^-) d- thus /(x) = ^ j e-^m ^e 

In the text, we systematically drop the constants such as ,3^2 since they are not relevant. 

• The Fourier multiplier with symbol m(^) is defined by 

m{D)f = F"^ [mFf] . 

• The bilinear pseudo-product with symbol m{^, rj) is given by its Fourier transform 

J" [Tmif, 9)] (0 = I MC, ri)?{il)g{i - ri) dv- 
Similarly, the trilinear pseudo-product with symbol m{S^, r], u) is given by 

^ [Tmif, g, h)] (^) = / m(^, 7], i')f{v)g{r])h{i - r] - u) dr] dv. 



H^ is given by the norm ||/||^^ = \\{D)^ f^. 
W^'P is given by the norm ||/||h/''.p = \\{D)'^f\\p. 



3 A formulation adapted to energy estimates 



Our aim here is to rewrite the equation in such a way that its dispersive properties become more 
transparent, but energy estimates can also be easily obtained. 
Spht 
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where Vp contains the unknowns which (in the linearization (jl.3p ) propagate as electromagnetic 
waves, and Va the unknowns which (still in the linearization ()1.3|) ) propagate as acoustic waves. 



3.1 The acoustic system 

We focus here on the evolution of Va = {Qu , n , QE ,0). It is governed by the system 

dtQu = -QE - V^ - c^pVp 
dtn = —V • {pu) 
V ■ E = -n. 

In order to diagonalize this system, let us switch to the unknown function 

^~ 2 \ \D\ ^^ \D\ ^ 



so that 



Qu 



V \B\ 

-2— -am^ and n = 2 ' ' IHe^. 



The evolution of A is given by 



2dtA = 2i{D),^A - ^^^ • {nu) + ^ f |u| V c^W' 
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3.2 The Maxweh (or electromagnetic) system 

We focus here on the evolution of V^ = {Pu , , PE , B). By (jl.2p . it suffices to consider PE and 
B. These fields are governed by the equations 

dtB = -V X E 

dtPE = V X 5 + P{pu) 

which implies 

d^B -AB + B = -Vx {nu). 

Setting 



it satisfies 
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and the original unknown functions Pu, PE and B can be recovered by 
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3.3 Summarizing 

The Euler-Maxwell system now reads 

' dtA - i{D),^A = -i W • {nu) + \i\D\ (\u\^ + |n|2 



{EM') 



dtB-i{D)B = -^^x{nu) 
[ iA,B){t = 0) = {Ao,Bo) 



with 



Qu 
n = 
Pu 



V r. 



, \D\ 



-m 



3mA 



mtA 






X 3mB. 



The data {Aq,Bq) of [EM') are easily expressed in terms of the data {uq^uq^Eq^Bq) of {EM): 
A =^^(^no + i^-^o) and B^ = -^x E^ + t^-^B^. 



Let us finahy define the profiles of A and B 



a{t) 



d^f ^~it(D), 



A{t) and b{t) =^ e 



^^f p-it{D) 
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4 A formulation adapted to decay estimates 

As we saw, the system {EM') written above is equivalent to {EM); it will be the correct formulation 
to perform energy estimates. However, as far as dispersive estimates go, we will not need all the 
structure of {EM'): only resonances will play an important role. It will be convenient to write 
{EM') in a more compact form. 



4.1 Duhamel's formula in Fourier space 

Writing Duhamel's formula in terms of a and b gives 



a{t) = Ao + Jq e 
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Taking the Fourier transform gives 



{a{t,^) = Ao{^) + "nonlinear term" 
b{t,0 =Bo{0 + "nonlinear term". 

In order to make notations lighter and estimates easier, we will now give up some of the structure 
of the above system. 



Convention 1. We will denote indifferently C{t) for A{t) or B{t), or their complex conjugates, 
and c{t) for a{t) or b{t), or their complex conjugates. Similarly, we denote e^**' " for any of the 
groups e^*<^>, e-^*<^>, e'^^^^'^' , or e-^*<^>'== . 

For instance, u or n are a linear combination of terms of the type e^**'^)*c. 

It is always understood that, in an expression of the form e^^*'^^''c{t), the meaning o/ e^**'^" 
is consistent with that of c. For instance, if c stands for a, then e^**'^/* stands for e**'^'''^ . 



With this convention, it is easy to see from the above that the "nonhnear terms" can all be 
written as a linear combination of terms of the following type (which we denote generically by g) 



git,C) 



5«s0(€,r;)^^^^ Oc(??, s)c{^ - r/, s) dij ds, 
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where m is such that 



dfmoiC) 
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"^(C)^) = '^o(C)™'i(^)'^2(C ~ ^) with < 
and (p is one of the 4>]^fl^ given by 

where 61,62 = ± and k,i,m are either 1 or c^. 

4.2 Space-time resonances in the context of Euler-Maxwell 

Seeing (j4.1|) as an oscillatory integral, it becomes clear that the cancellation properties of (j) and 
drf(j) will provide a key to understanding the large time behaviour of our system: this is the idea 
of space time resonances. See ^llj for a general presentation, and [lOj for the case of (semilinear) 
Klein-Gordon equations with different propagation speeds. 

Recall that the phase functions corresponding to all possible quadratic interactions are given 
by 

'PITJ^^ ^) '= iOk + ei(r?). + 62(^ - vU 
Next define for each interaction the space, time, and space-time resonant sets 

def 



(4.3) 






{{i^ V) I <Pl]i% = 0} ("space resonances") 
kim - { (?' ^) I (^r,(t)^kfm = ^i ("time resonances" ) 






7^' 



ei,e2 def ^ei,e2 ^ T-ei,<:2 
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The set of all space-time resonances is 



nr; 



("space-time resonances") 



']-> II ■T-)eiie2 . 

'^ ~ ^>^i,e2,k,e,m''^k,e,m^ 

it is compact and hence it is bounded. We denote by C-ji — 1 the radius of a ball that contains TZ. 
Finally, define the outcome and germ, or source frequencies for space-time resonances: these are 
simply the projections of TVj^f^ in the ^ variable, respectively the union of the projections in the r] 
and i — T] variables. More precisely if vr^(^', rj') = ^', 'Kr^H' ■, f]') = v' ^iid 7rg_^(.^', r]') = i' — rj' , we set 

o =^ TTg (7^) 

G =^ vr^ (7^) U 7rg_^ {IZ) . 

Definition 4.1. Space-time resonances are said to be separated if no outcome frequency is also a 
germ frequency. In mathematical terms, Q DO = 0. 
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5 Some linear and bilinear cutoff Fourier multipliers 

We overtake here some of the cut-off functions defined in |1U] : see Proposition 112. Il for results on 
the boundedness of the associated operators. 

5.1 Low or high frequency cutoff: Zi, Z^ 

First pick Mq large enough (the precise value of Mq will be fixed in the following, for the moment 
it is simply > C-ji defined in Section 14. 2|) . 

It will be necessary in the proof to distinguish between high and low frequencies. To this end, 
we introduce the cut off function 9{^,rj), which is such that 

eeC^ , = Ion 5(0,1) and = OonS(O,2)^ (5.1) 

Restricting to high, respectively low frequencies, is achieved by the operators 



5.2 Cutoff for O: the operators Zq, Z^ 



o 



Recall that O and Q were defined in Section | 

Under the resonance separation condition (definition 14. ip . it is possible to find 6q such that no 
frequency in BiQSf^{0) (a lO^o-neighbourhood of O) is a source of a space-time resonance: 

5lO5o(O)Ug = 0. 

Define xo a smooth cut-off function such that 

Xo = 1 on Bs^/2{0) 

Xo = outside of Bs^^ (O) 

Xo + Xo = 1- 



and then let xo satisfy 



The corresponding operators are 

Zo = Xo{D) and Zq = Xo{D). 

5.3 Cutoff for S and T: the symbols xs ^ind xt 

The cut-off functions which we are about to define will, for a given set of indices ei,e2,k,i,m 
separate 7^*^^'^ from S^^f^; of course this can only be done away from a neighbourhood of 7^^^/^, 
where these two sets intersect. Dropping for simplicity the indices, the function xs localizes away 
from T, whereas XT localizes away from S. Since T = {(j) = 0} whereas S = {drjcf) = 0}, this 
explains the inequalities ()5.2p . 



Lemma 5.1. For each set of indices ei,e2,k,i,m, it is possible to find cut-off functions 

Xjfi'^2{C,r]) , X5,^i/2 (C, ?/) 
such that (in the following, we drop the indices ei,€2,k,i,m for simplicity) 
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XT; XS 0.1^^ smooth. 

Their sum equals one away from IZ: XT + X5 = 1 ^/dist((^,r/),7^) > 5q/W. 



The derivatives of ^ and 
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if\a\ < 20, then 
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(5.2) 



for some integer uq . 



5.4 Paraproduct decomposition: the symbols (^ and ^^ 

Following the original idea of Bony [3j, we would like to distinguish between regions where \r]\ > 
\C-v\, respectively |^ - r/| > \r]\. 

We pick two functions C^{^,r]) and C^(^)^) such that 

• C^ and (^ are smooth. 

• C^ and C^ are homogeneous of degree zero outside of 5(0, 1). 

• C'^{^,v)+C\tv) = 1 for any (^,r/). 

• If KC,"^)! > 1 and (^,7]) € Supp(C"'^), then |^ — ??| < c\r]\ for a constant c. 

• If K'^,^)! > 1 and (^,7?) G Supp(C^), then |?7| < c\^ — r]\ for a constant c. 



6 The a priori estimates and plan of the proof 

In order to prove Theorem II. H we will prove the following a priori estimates, valid if e is small 
enough. 

Energy estimate 

• ||(^, ;B)||j|^jv < e{t) ""^ for a constant Co, and any t (regularity in L^). 

Decay estimates 

e 



\\iA,B)\\ „„,, 



w 



ii-'^y ~ (t)i+35i 



w 



ZoiAB) 

^oiAB) 
Localization estimates 



-1 < 



(i+^i) ~ (t)l-3<5i 



(square integrable decay above L ). 



(decay slightly below L^ for "non-outcome" frequencies). 



Zo {u, n] 



< -j-T (decay ~ j in L°° for "non-outcome" frequencies). 



|x|(a, 6)||j:^jv' < ^v(f) (localization in H^ ) 



\x\ ' Zc){a,h) < e (localization in L? for "non-outcome" frequencies). 
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The constants N, N' , N" are chosen such that N - Ni > N" - Ni > N' > Ni, for a sufficiently big 
constant A^i; in particular, N is sufficiently big for the necessary arguments in [TU] to apply. The 
constant 6i is chosen sufficiently small for the necessary parts of the argument in [10] to apply. We 
call II • ||x the norm corresponding to the above quantities: 



\\{A,B)\\x=^suj> 
t 



(()-''"■ l|{AB)ll„» + {*>*+"' ll(AS)|| 



„,»~,(^..)- +<*>""' «<^-'"«...» 



w 



".(^+^1 



+ (t) Zo{A,B,u,n) 



+ 



|x|(a,6)||^^/+ \x\^/^Zo{a,h) 



Since local well posedness is easily dealt with, and the data are chosen such that 

||(e^*<^>-A,e**<^>i3o)|U<e, 
the proof of the theorem consists in proving the following a priori estimate: 



\\{AB)\\x<W'^^^^^Ao,e^'^^^B,)\\x + \\{AB)\\\ + \\{AB)\\ 



X- 



We will proceed by showing that all the quantities appearing in the definition of X can be controlled 
by the above right-hand side. More precisely, the plan is as follows 

• Decay estimates are proved in Section [71 

• Localization estimates are proved in Section [HI 

• The energy estimate for A: sup^(t)~*-^'"^||^(t)||j:^jv < e is proved in Section [TOl 

• The energy estimate for B: sup^(t)~'^"^||;B(t)||j:^jv < e is proved in Section[9l 

• Finally, in section [11] we give a sketch of the proof of the scattering. 

7 Decay estimates 

We want to prove here that 



sup 

t 



{t)l+''^\\{A,B)\\ 



w 



[^-"h 



,-i+{t) 



l-35i 



ZoiA,B) 



w 



^",(i+^i)" 



i+{t) ZoiA,B,u,n] 



< 



3^*<^>-A + e'*<^>So) +\\iA,B)fx. 



(7.1) 



7.1 Control of the W 



N" 



\+Si 



" and W^"'(^-^0" 



norms 



The two first norms in (|7.ip above can be controlled as in [lUj : 

• As far as the control of the W ''-6+ ij norm goes, the main difference between the Euler- 
Maxwell system and the setting of [10] is the quasilinearity of Euler-Maxwell. This induces 
a further loss of regularity in the nonlinear term, which is however easily absorbed using the 

ttN 

H norm. 

• The estimate for the W 'V3~ iJ norm is a low frequency question (since it is only prob- 
lematic on O). Therefore, the argument of [lU] applies identically. 

We do not detail these two points, and focus directly on the third norm in (|7.1|) . 
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7.2 Control of the W^'°° norm 



Proceeding as in Subsection 14. 11 we can derive a generic term g corresponding to the nonlinear term 
in Duhamel's formula for u and n. It turns out, since u and n are given from A and B by the action 
of a Fourier multiplier, that this g would satisfy exactly the properties listed in Subsection 14. 11 
Thus all we need to do is to prove that, for g as in Subsection 14. H 



„,.„^(!>ll(AS)lli. 



Zoe"<''>'9(*) 
In order to prove this, we shall split Zqb^^^ '^g{t) as follows 

^(Zoe^*<^>^g(t)) = / / xo(0e'^'^m(e, v)c{s, v)cis, ? - r/) dr? ds 

XoiOxsi^, ??)e*'"^m(C, r])c{s, ■n)c{s, ^ - v) drj ds 



+ 



+ 



XoiOxri^, vV'^'mi^, r])c{s, r])c{s, ^-rj)d'n ds 



(7.2a) 
(7.2b) 
(7.2c) 



In the above, we have used the cut-off functions xs and XT- Remember that these were defined 
in 15.11 depending on the quadratic interaction considered; they were therefore labeled Xq*.',"! and 



-'S1.S2 



XT-fc,!,m- The above equation is written in generic terms, but it is tacitly understood that the cut-off 
functions used are the ones corresponding to the quadratic interaction at hand. 

7.3 Preliminary estimate on dgC 

Observe from subsection 14.11 that e*''^^''=9sc(^) can be written as a sum of terms of the type 



m{i,r])C{r])C{i - T]) ds, 



where m satisfies the estimates of that section. Therefore, by proposition 112. H 



< 111 ||2 



(7.3) 



7.4 The small time term ( I7.2al ) 



Using repeatedly the Sobolev embedding theorem, and the dispersive estimate ()12.3p gives (assum- 
ing t > 1, the case t < 1 being trivial) 



e''^^^'T-^^M 



W^'' 



5^(*-^)<^)^T^^(^)^(^,,)(C,C)ds 



W'--'''^ 



<1 



^xo(0™{C:'y)(*^'*^)'^'^ 



Vl/5.6/5 



ds 



~ 7 / ll'-llvy6,i2/5 

^ 1 /" \\r>\\2 r ^ 1 \\n\\2 

5-, - / L- H-7 as 5- - C k- 







t 
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7.5 The term away from T (I7.2bp 

In order to deal with this term, integrate by parts in time using the identity jzdge^'^'^ = e^^"^ ■ Thus 



1 



^M= I Xo{Oxs{^,v)m{tv)j^C{t,r])C{t,^ - 7j)dr] 

XoiOxsi^, r])m{^, r]) — e'^'t'dsc{s, r])c{s, ^-rj)d'n ds 
+ {symmetric and easier terms}, 



(7.4a) 

(7.4b) 
(7.4c) 



where the "symmetric and easier terms" correspond to the case where the partial derivative dg 
hits the other c, and to the boundary term at s = 1. Using successively the Sobolev embedding 
theorem 112.11 and Proposition 112.11 gives 



git(i3).^-i([73^ 



VK2.' 



< 






icj) 



w^ 



w^-y^-'"^) 



< \\c\\ \\c\\ 



~ ^iTe^rll'^ll^- 



In order to estimate ()7.4bp . split it as follows 
■F-^ (l7:4b]l = 



/ Xo(6x<s(?, v)ni{i, r]) — e''"^dsc{s, it])c{s, i - r]) dn] ds 
1 J ^4> 

t-l rt 



+ 1 ...'Ml + IL 

t-i 



Use the Sobolev embedding theorem, the dispersive estimate p2.3p . Proposition 112.11 and the 
preliminary estimate (|7.3p to get, for (5 > small enough. 



e^'^'^^^F-^I 



< 



< 



t-1 



r 



(t _ s)(3/2-35) 



T~ 



Xo(Oxs(i,v)m(^,V) 7^ 



1 e 



i^^'(^)'{dsc),c 



iy5'(i-'5)~^ 



ds 



1 (t - s)3/2- 



3<5 



e— "(a,c) 



t-i 



1 



1 1 



iy5+n,3/2 I|C|Ih/^+-,(1/3-^)-1 '^S 



^ (t_s)3/2-35 5 5l/2+3<5ll^ll^^^ 



^ ll*^ 11^^3/2-35 
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As for //, use repeatedly the Sobolev embedding theorem ll2.H Proposition ll2.1l and the prehminary 
estimate (|7.3|) to get 



^it{D),^-ljj- 



W 



< II t^^ttW 

2,cx, ~ Ir ■'^llvt'4,2 



< 



< 



< 



< 



t-i 
t 



t-i 
t 



'^xo(i)xs{i,v)Hi,v)j, ' ^ 






14^4,2 



ds 



Vl/6.(l-'5)" 



(is 



,±is(D>, 



(5.c) 



t-1 
* 1 1 



i;^A6+n,3/2 



||C||^6+n,(l/3-<5)-l ds 



f 1 1 11^112 - 

Jt-l SsV2+3<5ll 



< iiri|2 ^ 



7.6 The term away from S ( I7.2cp 

First transform this term by an integration by parts using the identity . iq iy2 " ^jjc'*'^ = e*'^'^. This 



gives 



UM 



a, 

is\drj(l)\ 

a, 

1 J is\drj(j)\ 

+ {symmetric and easier terms} 



XoiS.)xrit V)-J^-Ji2^^^' ri)e''"^dr^c{ri)c{i - rj) dr] ds 



XoiOxrit V)-J^-Ji2^v'^^^' ri)e''"^c{ri)c{i - rj) dr] ds 



(7.5a) 

(7.5b) 
(7.5c) 



Let us begin with ()7.5ap . which we split as follows: 



- (I73il) 



d: 



XoiOxri^., V) ■ |o .|2 "^(^' v)e''^dr^c{v)ci(. - v) dr] ds 



IS 



d, 



t/2 i-t 

+ 

1 Jt/2 



def 



(7.6) 



/ + //. 



For 6 < 0, \6\ small, apply successively the Sobolev embedding theorem, the dispersive esti- 
mate ()12.3p . and Proposition 112.11 to get 



— 1 



< 



I 11(4^2,00 r^ 



< 



t/2 



t/2 



S Xo{0XT(i,v)-7Tf^rn{^,v) \ ^ ' 



1 1 



1 (t - s)--- s 



< 



t/2 



1 1 



i|9,;0|- 



Vk'-(h+*)" 



VK 



^■(i-*) 



1 ds 
^ (is 



/j. M •^A ll-^*^ H""+5 K' I r / 1 i'l^l ds 

I {t-sy^-^^s " w 

t/2 



n + 5,{i^-Sy 



(7.7) 



^L (^^ip^iii^ii^^;!^^^ 



< IICI 



1 



x: 
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Taking this time 5 > and small, and retracing the above steps, one gets 



IT-^II 



I < iirii^ 



The term ()7.5bp can be estimated in a very similar way. Indeed, since m satisfies the estimates ()4.2p . 
dr^m{^, 7]) yields at worst singularities of the type -pr, tt^-t • The above scheme can then be employed 



since by Hardy's inequality, and Plancherel's equality, 

8 Localization estimates 



wi^iO 



< wd^mh 



\xc\ 



We want to prove here that 
1 , 



sup 

t 



Vt 



X (a, 



\H 



N' + 



\x\^/^Zoia,b) 



< 



(e**<^>-A,e^*<^>^o) +ll(AS)llx 



As above, this reduces to proving that the generic term 

Jo 
defined in (|4.ip satisfies the localization estimates 



sup 

t 



^WmgllHN' + 



x\^/^Zo9 



\u,m 



i.i) 



By symmetry, it suffices to control 

9'it)=T-' [ e-'^«'^)m(r/,e)C'(C,r?)c(r?,s)c(C-ry,s)ds 
Jo 

(where the cut-off symbol (^, defined in Section [SJ ensures that \C — "hI ^ \v\ ^^r (^,??) large). The 
bound for the second norm in ()8.ip was derived in [10], and the same scheme of proof applies here 
(once again, the novelty compared to [lOj is that the Euler-Maxwell system is quasilinear, but the 
resulting loss of regularity in the nonlinear term is easily absorbed by the H norm). Therefore, 
we focus on the first norm in (|8.ip . for which some new difficulties arise. It will be helpful to split 
m as m = momim2 (see Subsection 14. ip . Multiplying g by the weight x corresponds in Fourier 
space to differentiating g with respect to ^. This gives 



d^g'iO 



+ 



+ 



+ 







e*"^m(C, ri)C^ {^, r])c{r})d{c{i - r/) dr] ds 
isd^(j)e^^'^m{^, rfjC} {S,, ?/)c(??)c(.f — r?) dr] ds 
e^'^'^mo{S,)'mi{r])d^m2{£, - r])C^{S,,r])c{r])c{S, - rj) dij ds 
e^'^^mo{S,)mi{r])m2iS, - 'n)9iC^iS,,v)c{ri)ciC - v) dr] ds 



+ d^moiS,) I I e*'"^mi(r/)m2(C - v)C^i^,v)c{r])c{( - r/) dr]ds. 



B.2a) 
^.2b) 
3.2c) 
^.2d) 
3.2e) 



16 



8.1 Estimate of (K2ai^ 

To estimate (|8.2ap , use the Strichartz estimate (|12.4p and Proposition 112.11 to get 

ft 



LF^^dHSSD 



fjN' r^ 



^ e-<^>r^(5,,)^i(^,,)(C,e-<^>xc)ds 



HN' 



< 



Tmi^,n)(H^,r,){C,e''^^^xc) 



< 


l|C|| 




w 


< lic||^ 



i+iV' + l,(^-«l) 



-1 xc 



Ms.r^^'Mi-^.)-^ 



,(i+i^i)"' 



is) is) 



4~3<5i 



(i + i^i) 



< 



X 



(t). 



8.2 Estimate of (K2R\ 



To estimate ()8.2bp . distinguish three types of interactions, by writing c = Zqc + Zqc. The 
term (|8.2ep can be written as 



r 



-1/ 



2eD 



+ I e""^ 'sT^^^.^)^i(^^j^)Q^^{ZoC,ZoC)ds 



Js{D)^rp 



+ / e'°''^'sT^(^^,^)i^i^^,^)Q^^{ZoC,ZoC)ds. 



B.3a) 
^.3b) 
3.3c) 
^.3d) 



The term ()8.3ap can be treated exactly as in [lOj . thus we skip it. Next we shall bound the 
term (j8.3cp . The term (|8.3bp is comparatively easier, since the two interacting waves correspond 
to non-outcome frequencies, thus enjoying better bounds. As for the term ()8.3dp it is also easier: 
indeed for this term, the symbol C^iC^rj) imposes that ZqC is lower frequency than Z^C; but this 
is possible only if both are low frequency. 

Coming back to (|8.3cp . use Proposition 112.11 to get 



llMi3cj|^jv/ ^ s e*''^-^'r^(g^^)^i(g^^)a^0(Zc»C,Zc)C) 



HN' 



ds 



< / s||ZoC|| .1 .-1 
Jo d^-'^) 



ZoC 



w 



iV' + n,(l+ii) 



ds 



-'0 S2 



1 



1 



+3<5i sl-3<5i 



ds<. 



X 



Vt. 



8.3 Estimate of (1872^ 

By ()4.2p . d^m2{£, — rf) can be bounded by |g_ , . Bounding by Hardy's inequality \J_ \ c{^ — rj) in 
L^ by 9gc(^ — rj) in L^, the estimate for (|8.2ap can be easily adapted. 



17 



8.4 Estimate of (lOdb 

Since d^C^{^,r]) does not have a singularity, this term is easy and we skip it. 

8.5 Estimate of (187261) 

By ()4.2p . 5^mo(0 can be bounded by 1 for high frequencies, and ttt for small frequencies. Forgetting 
about high frequencies, which are easily dealt with, we need to bound 



in H . By Hardy's inequality, it suffices to bound 







F ^d^ j / e''"^d^mi{r])m2{^ - ri)C^{S,,r])c{r])c{^ - r]) dr] ds 



in H'^ . But expanding the £, derivative yields terms similar to (j8.2ap ()8.2bp ()8.2cp ()8.2dp . which 
we have just seen how to estimate. 

9 Energy estimates for the Maxwell part 

We shall prove in this section that 

Mh^ ^ \\{Ao,Bo)\\hn + \\{A,B)\\l + r 7^ \\{A,B)\\h^ ds. (9.1) 

Jo {^/ 

Together with (jlO.ip . this will imply that 

The following observation will be crucial: it follows from their definition that A and B control the 
physical unknowns u and n as follows: 

WQuWhn < \\A\\hn 

ll^^llHiv+1 < ll^ll^iv (9-2) 

9.1 Preliminary estimate: dgU 

It follows from {EM') that 

Therefore by the product estimates ()12.2p . 



/ \i 'IT^I/I |2,i |2 

(nu) + i\lJ\ \u\ + \n\ 



< {\\n\\HN + \\u\\hn) {\\n\\^ + \\u\\^) 

1 
< 



(9.3) 



,|+35i 
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9.2 Distinction between outcome and non-outcome frequencies 

Consider the integral equation satisfied by b: 

V 



b{t) =Bo+ I e 

'0 



i^D)^ X {nu)ds. 



Split n and u into Zqu + Zon, respectively Zqu + Zqu. This gives 



b{t) =13o+ f e-*^<^>-^ X (ZonZou) ds 



+ / e 

'0 
JO 



(^) 



{D) 



(9.4a) 
(9.4b) 
(9.4c) 



The term (j9.4ap is easily estimated: using standard product laws, and the Sobolev embedding 
theorem: 



wmmw 



JJN 



< 



< 



/* 

Jo 


-is{D) V 

(D) 


X [Zcm Zqu) ds 


IJN 






/■* r 


















/ 


Zon 




Zqu 


+ 


Znn 




Zqu 




Jo L 




UN 




L°° 




L°° 




H'^l 



ds 



ll(A^)|lx 



u-^,m 



{u,n)\\jjN ds 



is) 
^ I -^\\iA,B)\\HNds 



10 is) 

For the term ()9.4bp . we take advantage of the frequency localization of Zon and Zou to write, 
with the help of Bernstein's inequality, 

ft 



iidnaEDii 



HN 



-is(D) 



V 

(D) 



X (Zon Zqu) ds 



HN 



< / \\ZonZou\\ 



- 1^ 1(1-25)- 



ds 



^ / \\Zo'n.\\fi_^^-i\\Zc)uh^_^yi ds 






<\\{AB)\\j, 

10 \^i 



■ds<\\{AB)\\l. 



9.3 Interactions between outcome and non-outcome frequencies 

Thus we now take a closer look at (j9.4cp . wich reads 



([OcD = / e 
/o 



-is{D) 



V 

W) 



X {Zqu ZqPu) ds 



+ I e"''^''^X X {ZonZoQu)ds 
Jo W) 

+ f e-''^''^^^x{ZonZou)ds 
Jo W) 



(9.5a) 
(9.5b) 
(9.5c) 
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The first term, (j9.5ap . can be estimated with the help of the Strichartz estimate (|12.4p and the 
standard product law 112. 2t 



lldSSiDII 



JJN 



t vv _ 

e-*"<^>— - X {ZonZoPu)ds 
i-^} 



JJN 



< 



< 



Zqu ZqPu 

\\Zon 



^{h+¥^r\{--¥.+^u%-^^v 



3-°! 



ZqPu 



H. 



N+ 



1 + \\Zon\\j^N+i 



ZqPu 



-a-'-^y 



.(i+i*i)"' 



<\\iAB)rx it) 
<\\{Am 



-i-35i 



-(^+iM"' 



fx 



The estimates for the terms ()9.5bp and ()9.5cp are exactly the same, changing the roles of n and 
Qu. We will only treat ()9.5bp . The term ()9.5bp can be decomposed as 



i^l>h\\ 



,-is{D) 



V 



X {Zqu ZoZiQu) ds 



' (D) 

t y 

+ / e-''^^'^j^^x{ZonZhQu)ds 





{D) 



(9.6a) 
(9.6b) 



(recall that Zi and Z^ have been defined in Section [5l and that Mq has been chosen so big that 
ZhZo = Zh). The term ()9.6ap can be estimated exactly as ()9.4bp : thus we skip it. We are left with 
()9.6bp . Recall now that n and Qu can be written as 



n 



it) = 2-H^9^ee^*<^>=^a(t) and Qu{t) = -2-^ame^*<^>^=a(t). 

\P^)cs l"l 



This implies that the Fourier transform of (|9.6bp can be written as a sum of terms of the type 



■F (|9:6bl) 



„is0(5,r?) 



xoiv) 1 - X 



i-n 



10 J \ \ M 

where, for simplicity, we denote undistinctly a for a and a, (j) has the form 

Ht i) = (0 + ei(?7)c, + e2(e - V)c, with ei, es G {±1} 
and rh{^,r]) is a (matrix- valued) symbol satisfying the estimates 



m{(,,r])a{ri,s)a{^ — r],s) dr] ds (9.7) 



d?d^mirj,0 



< 



1 



\{^,rj)\\-W\- 



A crucial point will be that, on the support of Xoiv) (^ ~ ^ { m^ ) ) ' since Mq is chosen big enough, 
1^1 >> \rj\ ~ 1, and (j) satisfies the inequalities 



QaQl3 1 



< 



|^||a| + |/3| + l- 
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Integrating by parts in ([9 



using the identity -r^dgC^^'^ = e^^'^, and denoting 



gives 



/w(C,f?) 



def 



(7?)(l-g(^))mfaO 



T S9Mh =- I I e^"'^(«'''V(^, v)dsa{ri, s)a{C - V, s) drj ds 
J'<^i^'V)^^^^ ^)2(^^ s)9,a(e - ri, s) di] ds 
+ / e^*"^(«''^V(e,r/)a(?7,t)S(e-??,t)dryds 
— / fi{S,,r])a{ri,0)a{^ — ri,0) drjds. 



(9.8a) 
(9.8b) 
(9.8c) 
(9.8d) 



The only difficult term is ()9.8bp : thus we skip the other ones and estimate it with the help of 
Proposition 112.11 



LF-^ dObli 



JJN 



3^"<^>T^(Ae^*<^>9,a(s))ds 



HN 



<|*||T^(^,e^*<^)a,a(s) 



H^ 



ds 



< 



oo\\9sa{s)\\HN-i ds 



ds<\\{A,B)\\l. 



10 Energy estimates for the acoustic part 

We shall prove in this section that 



HN < \\{Ao,Bo)\\hm + \\{A,B)\\l + r 7^ \\{A,B)\\h^ ds. 

Jo w/ 



Together with ()9.ip , this will imply that 



||(A^)ll^iv<et^«^ 



(10.1) 



10.1 The equation (KT\\ 

First we rewrite the evolution equation ()3.ip satisfied by A. We will use the notation D 
We start by expanding the first nonlinear terms appearing in (|3.ip 
We start by taking N derivatives of (j3.ip . We get 



IDT- 
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d^d^A = i{D)c,d^A - -D ■ {d^nu + nd^u) 



i\D\ 



N 






(10.2) 



+ iu . l^id^'u + ic^J-^d^'n + < 



\D\ 
~2 



i{D)c,d^A - u ■ Vd^A + in{D)c^d^A + R^ 



where the rest terms Rf^ consist of quadratic lower order terms. In particular, we have 



/?f = D ■ Id'^nu) - d'^nu + nd'^ u 

+ ffl (a^(|7x|2) - lud^'u + c^a^dnp) - Iclnd^'n) 

2R^ = 2i?f - [b, u]d"n - [b-,n]d^u 

+ i[\D\,u]d^u + icl[\D\,n]d^n 



nN _ nN 

tin — rtn 



m- 



1 .d^n + iv^I^d^'Pu := R^ 



2\D\~ ' ' ' 2 
It is clear that we have the following estimate 

\\R''\\l^<\\V{AB)\\lA{AB)\\hn. 

Hence R^ does not loose derivatives, but for the part of {A,B) which is out-come we do not have 
enough decay to use Gronwall directly. 



10.2 Non resonant phase 

Due to the slow decay of the Zqu and Zqu, we have to use non resonant properties of the second 
and third terms on the right-hand side of ()3.ip . 



Lemma 10.1. There exist a positive number kq > and a constant Cq > such that for |^| > Cq 
and |r/| < C-ji, we have 






1 



€im^v) 






(10.3) 



for ei, £2 = ± and k,£ = 1, c^. 



Proof We will only consider the phase (/>c3',cs,cs=(0cs + {v)cs ~ {^ ~ 'n)cs since the other phases are 
easier. Furthermore, we only prove the estimate on -r, not its derivatives. We have 



{^-V)c 






1 



cm' 
1 



^ \rj\^ 

o ^-v _^ \v\ 



2cm' ^m' ' mi'^^^ie^ 
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Hence, we see that {£, — r])c^ — Cs\^\ — Cs\r]\ < jfr, from which we deduce that 

C 



> iOc. + {V)cs - Cs\i\ - CsH 



\i[ 



Hence, if Cq is big enough then, 0t',c„c. > ^1+E^lL^i^ > q. 



D 



10.3 Energy estimates 

The Sobolev estimates for the Maxwell part were performed using simply Strichartz estimates and 
integration by parts in time depending on the cases. Due to the further loss of a derivative, this 
method does not apply here. Instead we will perform an iterated energy estimate that we find 
interesting in its own right. 

Using that u and n are both real, we deduce that 



a, 



Hence 



2 



9^e / V.u\d"A\^ + i[{D)^^,n]d^Ad^A + R^d^A. 



I 5He I {V.u\d^A\^ + i[{D)^^,n]d^Ad"A + R^d^A) ds. 



(10.4) 



(10.5) 



We would like now to explain how to control the three terms on the right-hand side of (jlO.Sp . 
For the first term, we split u into the outcome and non-outcome parts u = Zqu + Zqu. The 
non-outcome part has enough decay to apply directly the Gronwall argument. Hence, we will only 
concentrate on the outcome part. We recall that the profile a{t) associated to A is defined by 

A{t) = e*<-^>=«*a(t). Also, we have 



Zou = Zo 



V 



e^t{D)j^_^-u(D)-^ y 



a_e^(^)-.*a). 



(10.6) 



We denote by e^**'^"c(t) = C{t) the divergence of any one of the four terms appearing in (jl0.6p . 
To control the first term in the right-hand side of pO.Sp . it is enough to rewrite it in Fourier space. 
Hence, it is enough to consider 



>'^(^'^)m(e, ri)c{s, r])d^a{s, ^ - r])d^a{s, i)dr] d^ ds 



(10.7) 



^ def 

where the phase (j) is given by (fii^^rj) = (^ — r])c, — {C)cs =t (??)«and m{^,r]) = xoiv)- Split 



m{i,r]) = ^(lln )"^(^'^) + 



1-^'A 



m{^,rj). The first term corresponds to low frequencies of 



d^A, which are easily estimated; thus, we shall consider in the following that 



m{tv) 



Mo 



xoiv) 
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From Lemma [10. 11 we know that (p is always bounded away from zero in the support of m. Hence, 
we can integrate by parts in time (using the identity -Kdge^'^'f' = e*'"'^) in pU.Zp and get 



idinsD 



+ 



m{^,ri)——c{s,'r])dsid^a{s,^ — ri)d'^a{s,£,)jdrid£,ds 
'm{S,,r])——idsc{s,ri)d^a{s,^ — r])d^a{s,^)]dridS,ds 

it<j> 

— m(e, v)c{t, ri)d^a{t, ^ - r])d^a{t, i)dr] d^ 



1 



-m(^, 7?)^(0, 7?)9^a(0, ^ - r/)5^a(0, ^drj d^ 
We rewrite the time derivative in (|lU.8a|) as 



dsd^a{^ - r])d^a{i) + d^a{i - T])dsd^a{0 
From (|1U.2|) . we deduce that 

Hence, (jlO.Sap can be expanded as 



(10.8a) 

(10.8b) 

(10.8c) 
(10.8d) 

(10.9) 
(10.10) 



pOHal) 



1 



m{i,r])-C{ri) 



u ■ Vd^A{^ - r])d^A{0 + d^Ai - V)u ■ Vd^AiO 



m{i,i)-C{vi) 



m{i,rj)-C{ri) 



n 



{D),^d^A{i - ri)d^A{0 - d^A{i - ri)n{D),^d^A{0 



R^iC - v)d^AO - d^A{^ - 7])R^{() 



drj d^ ds 



dr] d^ ds 
(10.11a) 

dr] d^ ds 
(10.11b) 

(10.11c) 



The difficulty in bounding (jlO.llap is that A appears with A^ + 1 derivatives. The main idea is to 
use some sort of energy estimate to perform an integration by part so that the extra derivative can 
be moved on a term with fewer derivatives. Keeping in mind that u is real-valued, we can rewrite 
(llO.llap as 



(|l().llap 




mi 



(e, r])-C{ri) uiu) -{i-ri- u)d^ A{i - r? - v)d^ A{-e) 



+d^A{i - ri)u{u) ■ (-e - i^)d^Ai-C 



dr] d^ dv ds 




C{r])u{i^) ■ ]i{C, r]){i-r]- v)d^ A{i - r] - i^)d^Ai-0 

+/x(^ - u, r])d^{C -V- '^){-0d^{-0] dr] d^ dv ds 

where /i(^, r/) = "^(c) and we made the change of variable ^ — )• ^ — z^ in the last line. The integrand 
of the term appearing in the last two lines can be rewritten as 

C{r])u{u)d^{^ -T]- u)d^{-0 ■ []i[i, r])[i-r]-u) + ^(^ - u, r?)C] 
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and the term between brackets is equal to 

Mi^, V, v) = -K^, v)iv + '^) + Ht v) - K^ - z^, v)] ^■ 

Proposition 112.31 gives the desired conclusion, namely that 



(10.12) 



riO.lOa) 



l2^ 



< 



M{C, ri, u)C{ri)u{ij)d^A{^ - r? - u)d^A{-Odri d^ du 



ds 



\v\\ liril ll/9^>lll (is 

I "'II OO 11^ II OO r^ 9 



(10.13) 



<\\{AB)\\], 



1 



{s) 



1+65; 



-ds<\\{AB)\\\ 



The treatement of (|10.11bp is very similiar and we do not detail it here. 

To control the second term on the right-hand side of (jlO.Sp . we rewrite it in Fourier space. 
Hence, it is enough to consider 



2"*c(r?)((e - 7?),, - (0.ja^a(e - v)d''a{Odr^dids 



(10.14) 



where the phase 4> is given by (l){^,rj) = (^ — 1])^^ — {^)c^ ± {ri)c^ and e^^^^'^^'^c is one of the two 
terms appearing in the decomposition of n as n = e'^^^''">^N + e~'^^^''''>^N . The estimate of (jl0.14p 
is exactly the same as the estimate of ()10.7p and we do not detail it again. 

Now, it remains to control the last term on the right-hand side of (jlO.Sp . namely the term 
involving the rest term. Again, if the low frequency term is non-outcome then we can estimate 
the term directly using the integrable L® decay of the non-outcome part. Hence, the only difficult 
terms are those for which the low frequency term is outcome. The most difficult terms are very 
similar to those we treated above by integration by parts in time, namely ()10.7p and ()10.14p . In 
addition we have terms of the type 



e'"^c{'q)d^ a{i - T])d^a{C)dr]d^ds 



(10.15) 



where the phase (p is given by 0(^, rj) = —{S,— r])c^ — {(,)cs i {'r])cs and terms obtained by taking the 
complex conjugate These two types of terms are even better than ()10.7p since integration by parts 
in time gains a factor t4t that behaves like ^ in the dangerous region |r/| << |^|. 
We also have terms of the form 

ft 



^'^c(r?)5^-'=a(e - r])d^ a{i)diidids 



(10.16) 



where A; > 1 and the phase (j) is given by 4>{S,,ri) = ^{^ ~ 'n)cs =t {S,)cs =t (^)« a-iicl e^**^-^/'c(t) = C{t) 
denote outcome (low frequency) terms and we denote undistinctly a for a or a or their complex 
conjugate. From Lemma llO.ll we know that (p is always bounded away from zero in the region we 
are interested in, nemaly ^ large and |?7| < Cti. Integration by parts in time yields terms that are 
easier to control than above. In particular the corresponding term to (jlO.Sap can be expressed as 



^C{v) 






u ■ V5^-M(e - v)d^AO + 5^-M(C - v)u ■ Vd^AiO 



n{D)^^d^-kA{i - ri)d^A{i) - 5^-M(e - r])n{D)^^dNA{0 



RN-k^^ _ rj)d^A{0 - d^A^ - v)R^{0 



(10.17) 
(10.18) 
(10.19) 
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which can be easily estimated. 

Finahy, we also have terms for which Pu carries the greatest number of derivatives. For these 
terms, we cannot use the cancellation coming from the energy estimate. To gain the two factors of 
1^1, we take advantage of the fact that Pu is more regular, namely it is in H^~^^ and the fact that 
the phase (p involved in this case is bounded below by |^|/C. The term corresponding to ()10.7p is 
of the form 



e'''^c{ri)d^+^{^ - ri)d^a{Odvd^ds (10.20) 

10 J J 

where the phase 4> is given by 4>{S,,'r]) = i(^ — ^)i — iOcs ^ iv)! ^^'^ e^**^^''c(t) = C{t) denote 
outcome (low frequency) terms. It is clear that \(p\ > \S,\/C in the region we are interested in, 
nemaly ^ large and |ry| < C-ji. Hence, we can perform an integration by parts in time and conclude 
as before. 

11 Scattering 

Let us prove for instance that A scatters. We write symbolically the equation ()3.ip on A as 

dtA-i{D)^^A = dCC 
By definition, A will scatter in ff , say at +oo, if and only if 

/ e''^^'^^^dC{s)C{s)ds 
Jo 

converges as i — )• oo. By the Strichartz estimates ()12.4p . it suffices that the right-hand side dCC 
belongs to L^^ ^ i[0,Qo),Lx ]■ This is the case since 



"'"'"" (^^iM-\(4-o-^ -"'""- ^^"'^y^-''' .(i.!M-^ <^' 



where the last inequality follows since e is small enough. 



12 Appendix: analytical tools 

12.1 Sobolev embedding theorem 

lil<p<q<oo and 



then 

12.2 Product laws 

If 1 < p, r < oo, 1 < q < c«, k > and 



7 3 3 

P Q 



q ^ \\J WW^'l 



1 1 _ 1 

p q r ' 



(12.1) 



then 

WfdWw-r^WfWw-rMg + WfWJgWw'^.r (12.2) 
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12.3 Dispersive and Strichartz estimates 

The standard dispersive estimates for Klein-Gordon can be found in Ginibre and Vefo [T5 



MD) 



<t^ ^11/11 .(11) W if2<p<ooande>0. 



(12.3) 



We need the fohowing Strichartz estimate for the Klein-Gordon equation (see for instance Ibrahim, 



Masmoudi and Nakanishi [20]): if e > and < 5 < |, 



e'''^^^F{s)ds 







< II-FI 



iw^^Vw[i-¥+^)ii-^r 



(12.4) 



2 L\^^^") w\ 

For the reader familiar with Besov spaces, this estimate follows from the interpolation between 



e''^^'^F(s)ds 







< IIFI 



L1L2 



and 



e'''^^^F{s)ds 



< \\F\ 



r2 d5/6 
^ -^6/5, 2 



12.4 Boundedness of multilinear Fourier multipliers 

After cutting off with the help of the functions defined in Section [5l the manipulations which we 
perform lead to various pseudo product operators. Their boundedness properties are stated in the 
following proposition. Notice that the statement below is very far from optimal, but sufficient for 
our purposes. 



Proposition 12.1. Assume that m satisfies the estimates (fj^. 

1 _ 1 I 1 



(i) Then for any p, q, r in (1, oo) satisfying - = — | — , and A; > 0, 



|2^m(/)5)llH'fe,'- ~ ll/llvK'=+i'P IIS'lliy'^+1'9 



(a) Assume 



1 ^ J\ 

Then there exists a constant, which we denote n > 0, such that for any p,q,r in (l,oo) satisfying 
1 = 1 + 1, andk>0, 



\TM 



p-\j jyjww'' 



< 



y[/k + n,p \\g\\\Yk + n,q 



(Hi) Assume 



^(?, V) = "!(?, vX^it i) or m(^, r/)C^(^, r?)5g</>(^, r/). 



Then there exists a constant, which we stil denote n > 0, such that for any p,q,r in (1, oo) satisfying 



i = i + i, andk> 0, 



(iv) Assume 



/i(^, r/) = m{^, r/)C^(^, r]) or m{^, r/)C^(^, r])d^4>{^, r/). 
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Then there exists a constant, which we stil denote n > 0, such that for anyp,q,r in (l,oo) satisfying 
1 = 1 + 1, andk>0, 



< 



(v) Assume 
where m and -r satisfy the estimates 



^/i(/;fl')llvyft,'- ^ ll/llp \\9\\w''+"'P ■ 



H^,v) 



|9f5,^m(e,7?)|<(|e| + |r?|)-H-l/^l and 



Then 



dp^- 



\T^^if,g)\\Hk < ll/lloolblbfe+i 



<m + \v\y 



H-m-1 



Proof. Estimates similar to the first four points above were proved in |10j . It essentially suffices to 
use the basic estimate 



\T,{1 



f-tyj 1 cj J w 'f 1—^ 



\f^\\H'^/2+e\\J\\p\\g\\q 



if e > and - + - = -, the estimates given in Section [S] on the various symbols, and a paraproduct 
decomposition to handle large frequencies. 

The fifth point follows from the classical Coifman- Meyer theorem [7J. D 

Next, we want to study a particular kind of symbol, which will not satisfy standard Coifman- 
Meyer bounds, but still admit Holder-like bounds (in the bilinear case for instance, we only focus on 
the case L°° x L^ ^- L^ bound, but it should be clear from the proof that more general LPxL^ — )■ U' 
bounds, with \ + \ = j.-, also hold). 

Lemma 12.2. Let R he a fixed constant 

(i) Let fi{^, r]) be a smooth symbol such that 

Supp/u C{\ri\<R} and d?d^f^{^,ri) 



<J_ 



for any ^, ij. 



Then\\T^{f,g)\\^<\\f\U\g\\2. 

(a) Let fi{^, r], v) he a smooth symbol such that 

Supp//C{|r?|<i?,|H<^|e|} and \d^dP^dZii{^,J],u 



< 



|^|l"l + l7l 



for any (,, rj, v. 



Then \\T^{f,g,h)\\^ < ||/||oo|blloo||/i||2- 

Proof. We take for simplicity R = 1, and first define standard Fourier space decompositions 

• Let C be a non-negative function, equal to 1 on i?(0, .9), zero outside of -6(0,2), and such 
that Ylj&s C(C - j) = 1 for any ^. Denote 

i-3 
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Let V' be a non-negative function, equal to 1 on -6(1, 1.5), zero outside of -B(.5,4), and such 
that J2jez3 V' ( ^ ) = 1 for any ^ / 0. Further denote 

m = E V' (I) and xio = E v^ (I) 

and the associated Fourier multipliers 



Proof of (i) Split n as follows 



The symbols ^j are uniformly controlled in C for any k. Thus they define operators which are 
uniformly bounded L°° x L^ — t- L^. Observe furthermore that, due to frequency localization 
properties, T^.{f,g) = T^.{f,Qjg); and that for the same reason, the families {T^.{f,g))j and 
{Qjg)j are almost orthogonal in L^. These arguments lead to the following inequalities 



\Uf^9)\\l < E \\TH(f^9)\\l < E 11/11- \\Qj9\\l < ll/IlL ll^lli, 



|2 ^ Y^ ll^l|2 11^ „||2 ^ ||^||2 II ||2 

loo I|y|l2! 

3 3 



provmg (i). 

Proof of [a) We will essentially run the original argument of Coifman and Meyer [7|. First set 
/i((^, r/, z/) = Jjl{£^ — rj — v,ri,v), and observe that the bounds on ^ translate into 

dtd',dZm,V,^)\<,^^y (12.5) 

Next split Jl as follows 

K^, ^' ^) = E "^ ( 2i ) ^'■^' ^' ^) = E /^J (^' '^' ^) ^^P *° ^ remainder, 
i i 

Since the remainder is compactly supported, and /i is closed, it will be easy to estimate, thus we 
forget about it and focus on the sum over j. The support of JijiCiVi^) is contained in a box 
{\v\ ^ 1 1 1^1 < '^''^^ ) I'^l ^ 2m2''^^}- I* ^an be expanded in (periodic) Fourier series adapted to the 
larger box {\r]\ < 2 , |^| < 2^^'^ , |z/| < 2^^'^ , and then recovered by restriction. This gives 

or, coming back to /i. 



/^i(^, ry, i^) = V' ( ^-^ j X(2r?)x ( ^j^ ) E "fc, 
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Next notice that the a'j^ ^ are uniformly bounded in j, with arbitrarily quickly decaying (inverse) 
polynomial bounds: 

for any N, sup la^^J = ak/,m < ^^j^j^^^^n ' (^^.G) 

This can be seen by simply coming back to their definition: 

y|^|<2J+2 JH<lJ\u\<2i-'^ 

J\^\<i J\'n\<i Jw\<i 

and the conclusion follows since the bounds (|12.5p imply a uniform control (in j) of the symbols 

/I,(2^+2^,,?,2-'-M- 

Coming back to physical space, we have achieved the following decomposition for T^: 

j,k,e,m 



The desired estimates follows easily by almost orthogonality between the j-summands 



\T,if,9,h)\\,< Y, 



k,i,m 



E 



a 



k,i,m 



< 



E 

k,£,m 



a- 



kl.mWJ Woo 



1/2 



1/2 



< 



E"M,-[ll/llLNlLl|/^lli 



211/2 



k,l,m 



< 



Iffllooll^lb, 



where we used in the last inequality the bound (jl2.6p . 

Equipped with the previous lemma, we can prove the following proposition. 
Proposition 12.3. Let M be as in \10.12\) . and fix a > 0. Then the following estimate holds: 

\\TM{f,9,h)\\2 < ||/||h/i+-,°°I|5||oo||/i||2- 



D 



Proof. Recall that 



def 



M(e,r/,z.) = -^(e,r/)(r/ + i/) + [^(C,7?)-^(e-i/,r?)]e = Mi(C,r?,z.)+M2(e,r/,z.) 



where the operator fi is supported in a strip {|ry| < 1 , |^| >> 1} and satisfies the bounds 



<J_ 



We will treat separately the operators Tmi and Tm2, further distinguishing for the latter between 
the regions where \£, — r] ~ '^\ ^ ^) ^^^ those where \^ — r] ~ ^\ >> 1^1) by writing 

Tm2 (/, g,h) = Y Tm2 {Pjf, 9, Sj+ioh) + ^ T^./^ (S'j-io/, g, Pj h) up to a remainder. 

j>0 j>0 
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Since the remainder is smooth and compactly supported, it is easily estimated, and we forget 
about it in the following in order to concentrate on the sum over j. Notice that we overtook the 
Littlewood-Paley operators Pj and Sj defined in Lemma 112.21 

The operator Tj\/^ . Simply observe that 

TMAf,9,h)=T^{Vf,gh)+T^{f,{Vg)h). 

Thus Lemma 112.21 gives the conclusion. 

The operator T^/j in the case |C ~ ^ ~ ^1 ^ ^- Recall that M2 is given by [/i(^, r/) — ^(^ — i/, r])] ^. 
There is no cancellation between the two summands in the range we consider, so /x(^, r/)^ and 
m(C ~ ^) 1])^ can be considered separately. Since they are estimated in similar ways, we focus on the 
first one. Notice that r„(f „)c = VT^. Using Lemma ll2.2l and proceeding in a straightforward way, 
we get the desired estimates: 



^VT^iPjf,gS,+ioh) 



<^2^ \\T,iP,f,gS,+wh)\\, <^2^\\Pjf\\oo\\9\\o.\\S,hh 
j j 

< ||/||;yl + a,cx) 11511211/1112- 



The operator Tm2 in the case \C ~ V ~ '^l >> ^- ^^ this case, we observe that the operator 

{f,9,h) ^ ^TM2{Sj-iof,g,Pjh) 



has a symbol 



M^(g,r?,z.)=M2(g,r?,z.)^V( ^ I " ] X ^ "" 



2J 



2J-10 



which can be written 

M^(e,r?,z.)=M2(e,r/,i^)-z. with M2(e,r/, i/) =^ J^ V' ( 



^ '^ ""^xi^Ml d^ii{i-iv,ri)dt. 



2i-io 
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The key observation is that, due to the hypotheses on //, we have the bound |C|'"'|9?)u| < C and 
hence M2 satisfies the conditions of Lemma 112.21 The estimate follows easily: 



||TM,(5,/,<7,Pi/i)||2= r^,(V/,<7,/i) <||V/|L||<7 



loo II"'II2' 
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